Effect of the 5-meson on the instabilities of nuclear matter under strong magnetic 

fields 
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We study the influence of the isovector-scalar meson on the spinodal instabilities and the distilla- 
tion effect in asymmetric non- homogenous nuclear matter under strong magnetic fields, of the order 
of 10 18 — 10 19 G. Relativistic nuclear models both with constant couplings (NLW) and with density 
dependent parameters (DDRH) are considered. A strong magnetic field can have large effects on 
the instability regions giving rise to bands of instability and wider unstable regions. It is shown that 
for neutron rich matter the inclusion of the 8 meson increases the size of the instability region for 
NLW models and decreases it for the DDRH models. The effect of the 8 meson on the transition 
density to homogeneous /3-equilibrium matter is discussed. The DDRH5 model predicts the smallest 
transition pressures, about half the values obtained for NL5. 

PACS numbers: 21.65.-f 26.60.Kp 26.60.-c 97.60.Jd 



Neutron stars with very strong magnetic fields of the 
order of 10 14 — 10 15 G are known as magnetars [1-3] and 
they are believed to be the sources of the intense gamma 
and X rays detected in 1979 (for a review refer to [4]). 
To date, 16 magnetars have been identified as short 7-ray 
repeaters or anomalous X-ray pulsars, although some are 
still unconfirmed candidates [5]. However, according to 
Ref. [6] , a fraction as high as 10% of the neutron star pop- 
ulation could be magnetars. These neutron stars are hot, 
young stars, ~ 1 kyear old. In [7], Thompson and Dun- 
can have considered turbulent dynamo amplifications in 
young neutron stars as a mechanism for generating the 
strong magnetic fields. 

In magnetars the crust is stressed by very strong forces 
which deform the crust and may crack it. Once the sur- 
face cracks, the violent motions blast particles along the 
magnetic fields, triggering gamma rays and x-rays. In 
particular, scientists believe that the giant burst of en- 
ergy experienced by the magnetar SGR 1806-20 in 2004 
was triggered by a " starquake" in the neutron star's crust 
that caused a catastrophic disruption in the magnetar 's 
magnetic field. SGR 1806-20 is the most magnetic ob- 
ject observed and has on the surface a magnetic field of 
intensity over 10 15 G [5]. 

In [8] the authors have shown how to obtain both the 
moment of inertia of a neutron star as well as the frac- 
tional moment of inertia of its crust only in terms of the 
mass and the radius of the star. For the fractional mo- 
ment of inertia of the crust has an additional dependence 
of the equation of state (EOS) which enters through the 
values of the pressure and the density at the crust core 
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transition [8, 9]. The transition density enters just as a 
correction to the fractional moment of inertia, while the 
transition pressure is the main EOS parameter which de- 
fines that quantity. Using this result, Link et al. [9] have 
obtained a lower limit for the neutron star radius with a 
given mass, from the glitches occurring in the Vela pulsar 
and in other six pulsars. This constraint will put severe 
restrictions on the acceptable equations of state of stel- 
lar matter if the radius and mass of a neutron star is 
measured. 

Spin-up glitches have also been observed in all known 
persistent AXP [10, 11]. However, different from the con- 
ventional low-field radio pulsar, the glitches in AXPs are 
accompanied by a significant recovery of the spin down 
rate of the pulsar. It would be interesting if information 
of the crust and star properties could also be obtained 
from these glitches. 

It was recently shown in Ref. [12] that a strong mag- 
netic field, of the order of 10 18 - 10 19 G, has large effects 
on the instability regions of nuclear matter. Relativistic 
nuclear models both with constant couplings and with 
density dependent parameters were considered. It was 
shown that a strong magnetic field can have large effects 
on the thermodynamic spinodal instabilities zones giving 
rise to bands of instability and wider unstable regions. 
As a consequence, it was predicted larger transition den- 
sities at the inner edge of the crust of compact stars with 
a strong magnetic field. The direction of instability gives 
rise to a very strong distillation effect if protons occupy 
only partially a Landau level. However, for almost full 
Landau levels an anti-distillation effect may occur. 

In this paper, which completes Ref. [12], we study the 
influence of the isovector scalar meson on the low density 
instabilities of asymmetric nuclear matter under strong 
magnetic fields. We estimate the density and pressure at 
the transition from the non-homogeneous to the homo- 
geneous phase in stellar matter from the crossing of the 



2 



EOS in /3-equilibrium matter with the spinodal and dis- 
cuss the effect of the magnetic field on the direction of in- 
stability. We consider two kinds of relativistic mean-field 
approaches: non-linear Walecka models (NLW) models 
with constant coupling parameters and density depen- 
dent relativistic hadronic (DDRH) models with density- 
dependent coupling parameters. The last models seem to 
give more realistic results at subsaturation densities [13]. 
The inclusion of the <5-meson brings to the isovector chan- 
nel the same symmetry existing already in the isoscalar 
channel with the meson pair [a, u>) responsible for sat- 
uration in RMF models [14, 15]. The presence of the 
(5-meson softens the symmetry energy at subsaturation 
densities and hardens it above saturation density, giving 
rise to stable compact stars with larger masses [16]. In 
[13] it was shown that the instability region in the isovec- 
tor direction becomes smaller if the 6 meson is included. 
Therefore, it is expected that in the presence of a strong 
magnetic field these differences become larger. 

In the present paper we consider the NLW models 
NLp and NL6 [15] and the DDRH models TW [17] and 
DDRH<5 [18]. Only NL<5 and DDRH6 include the 6- 
meson. 

For the description of the neutron star matter, we em- 
ploy the standard mean-field theory (MFT) approach. A 
complete set of the equations and the description of the 
method can be found in Ref. [12, 19]. The Lagrangian 
density of TW [17, 20] and DDRH<5 [18, 21] models reads: 
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where are the baryon (b=n, p) Dirac fields; a, to, 
p, and 5 represent the scalar, vector, isovector-vector 
and isovector-scalar meson fields, which are exchanged 
for the description of nuclear interactions and A 11 = 
(0, 0, Bx, 0) refers to a constant external magnetic field 
along the z-axis. The nucleon mass and isospin projec- 
tion for the protons and neutrons are denoted by m b 
and T3fc = ±1, respectively. The mesonic and electro- 
magnetic field strength tensors are given by their usual 
expressions: fl^ = d^uj u - 8 v lu^, P^ v = d^pu - d u p^ 
and F^ v = d^A v — d v A^. The nucleon anomalous mag- 
netic moments (AMM) are introduced via the coupling 
of the baryons to the electromagnetic field tensor with 
<V- = |[7/*)7w] an d strength n b with n n = g n /2 = 
-1.91315 for the neutron and k v = (g p /2 - 1) = 1.79285 
for the proton, respectively. The electromagnetic field 
is assumed to be externally generated (and thus has no 
associated field equation), and only frozen-field configu- 
rations will be considered. The density dependent strong 



100 



80 



1 I 1 I 

— DDRH8 
--■ TW 

■ NL5 
NLp 



T 



I 1 I 




I i I i I i I i I i I 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
P h ( fm " 3 ) 



FIG. 1: Symmetry energy for all models used in the present 
work. 



interaction couplings are denoted by T, the electromag- 
netic couplings by q and the nucleon, mesons masses by 
to. The density dependent coupling parameters are ad- 
justed in order to reproduce some of the nuclear matter 
bulk properties using the following parametrization 

Fi(p) =Ti(p sat )fi(x), i = a,uj,p,5 (2) 

where x = pj p sa u with 

l + h (x + drf . 

J l (x) = a l —2, i = a,w, (3) 

l + Ci (x + di) 



for TW, 



and, 



f p (x)=exp[-a p (x-l)], (4) 



fi(x) = a,exp [— bi{x — 1)] — c, (x — di) , i = p,8 (5) 

for DDRtW, with the values of the parameters mj, Ti, a,i, 
bi, Ci and di, i = a,u>, p,S given in Table I. 

In the sequel, we define the magnetic field in units of 
the electron critical field B° B = 4.414 x 10 13 G, so that 
B = B* B%. 

For the DDRtW model we use the parametrization 
given in [18, 21] except for the parameter r p (p sat ) which 
we increase so that the symmetry energy is 31 MeV at 
saturation density and not 25 MeV as in [18], see Ta- 
ble I. The inclusion of the S meson reduces the symme- 
try energy at subsaturation densities but makes it much 
stiffer at supra-saturation densities because the 6 meson 
field reaches saturation and the p-meson field increases 
always with density [14, 15]. This is seen in Fig. 1 where 
the symmetry energy is given for the all models we are 
studying. 

For the NLW models, i.e. NLp and NL5, we add to 
the Langrangian density, Eq. (1), with gi — Ti, the scalar 
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meson self-interaction terms 



^bm n (g a <r) 3 - ^c(g a af, 



where b and c are two dimensionless parameters. The 
coupling parameters are constant and given in Ref. [15]. 

The stability conditions for asymmetric nuclear mat- 
ter, keeping volume and temperature constant, are ob- 
tained from the free energy density T , imposing that this 
function is a convex function of the densities p v and p n , 
i.e. the symmetric matrix with the elements [22, 23] 



( d 2 T 
\dpidpj 



(6) 



is positive. At zero temperature, the free energy density 
coincides with the energy density. We define the thermo- 
dynamic spinodal at T = as the curve on the {p n ,P P ) 
plane for which the determinant of the Tij is zero. The 
eigenvalues of the stability matrix are given by 

A± = \ (Tr(jc-) ± v/Tr(^) 2 - 4Det(^)) 
and the eigenvectors 5p± by 
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In Fig. 2 we show the spinodal sections for several 
magnetic field intensities and for two NLW models: NLp 
without (5-meson (thin lines) and ~NLS including 5 (thick 
lines). As discussed in [12], due to the Landau quantiza- 
tion, the magnetic field has a strong effect not only on the 
size, but also on the shape of the spinodal zone. Unlike 
the B = case, the instability zone is no longer symmet- 
ric with respect to the p n = p P line. For the proton-rich 
matter, the inclusion of the (5-meson has only a small ef- 
fect on the spinodal, namely a small reduction similar 
to the effect already described for the B = [24]. For 
neutron-rich matter in the presence of a strong magnetic 
field, the effect of 5-meson is stronger: the NL5 model has 
a spinodal region larger than NLp. This effect is similar 
whether the AMM is included or not. 

The DDRH models behave in a different way due to 
density dependence of the coupling parameters. In Fig. 3, 



TABLE II: Predicted density, proton fraction and pressure at 
the inner edge of the crust of a compact star at zero temper- 
ature, as defined by the crossing between the thermodynamic 
instability region of npmatter and the /3-equilibrium EOS 
for homogeneous, neutrino-free stellar matter in the (p P ,p n ) 
plane. The AMM is not included. 



TABLE I: Parameters for the TW and DDRH5 models. These 
two models have the same parametrization for the a and to 
mesons. 



B* 


Model 


p« oss (fm- 3 ) 


y P 


P m (MeVfnr 3 ) 





NLp 


0.067 


0.013 


0.255 




NL5 


0.063 


0.010 


0.175 




TW 


0.085 


0.037 


0.523 




DDRH5 


0.086 


0.036 


0.269 


3 x 10 4 


NLp 


0.071 


0.043 


0.284 




NL<5 


0.071 


0.040 


0.264 




TW 


0.071 


0.056 


0.342 




DDRH5 


0.077 


0.056 


0.259 


5 x 10 4 


NLp 


0.084 


0.068 


0.467 




NL<5 


0.084 


0.064 


0.467 




TW 


0.082 


0.086 


0.461 




DDRH<5 


0.084 


0.086 


0.306 


10 5 


NLp 


0.105 


0.121 


0.822 




NL<5 


0.105 


0.117 


0.896 




TW 


0.101 


0.146 


0.673 




DDRH<5 


0.096 


0.148 


0.404 


2 x 10 5 


NLp 


0.128 


0.207 


1.203 




NL<5 


0.130 


0.203 


1.447 




TW 


0.128 


0.236 


1.016 




DDRH<5 


0.118 


0.237 


0.644 


5 x 10 5 


NLp 


0.162 


0.367 


1.275 




NL5 


0.170 


0.369 


1.830 




TW 


0.198 


0.402 


2.331 




DDRH(5 


0.228 


0.422 


4.460 



we present the results of the thermodynamic spinodal in- 
stability region for TW (thin lines) and DDRFM (thick 
lines). For B* = 10 5 and 2 x 10 5 the spinodal has three 
and two bands, see Fig. 3, corresponding to the occupa- 
tion of the first three and two Landau levels, respectively. 
For proton-rich matter with B* = 10 5 and 2 x 10 5 , the 
inclusion of the (5-meson reduces the size of the spinodal 
zone corresponding to the last occupied Landau level 
(LL), i.e., third LL for B* — 10 5 and second LL for 
B* = 2 x 10 5 . The other instabilities zones are not af- 
fected. For the neutron-rich matter, the (5-meson reduces 
the size of the spinodal zones for all the magnetic fields 
considered, contrary to what happens with NL<5. 

For B* = 5 x 10 5 the protons are totally spin polarized 
in all the models under study and for the complete insta- 
bility region. The size of the spinodal zone is the largest 
one. 

The existing differences between the different models 
for neutron rich matter will affect the transition density 
from a non-homogenous phase to a homogeneous phase 
in stellar matter. In fact, the density at the crossing of 
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FIG. 2: (Color online) Spinodal section in terms of p p versus p n for NLp (thin lines) and NL<5 (thick lines) at T = OMeV and 
for several values of magnetic fields (a) without and (b) with AMM. 
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FIG. 3: (Color online) Spinodal section in terms of p p versus p n for TW (thin lines) and DDRH5 (thick lines) models at 
T = OMeV and for several values of magnetic fields (a) without and (b) with AMM. 



the EoS for /^-equilibrium stellar matter with the ther- 
modynamic spinodal gives a reasonable prediction of the 
transition density [25, 26] of the crust to an homoge- 
neous phase in stellar matter. It was shown in [9] that 
the transition pressure and, at a second level, the tran- 
sition density define the fraction of the star's moment of 
inertia contained in the solid crust. In the sequel we will 
estimate for the different models and magnetic fields the 
transition density and transition pressure at the inner 
edge of a compact star in /3-equilibrium. We will discuss 
only the cases for which the crossing is occurring at the 
first LL. For fields B* < 3 x 10 4 at densities above the 
crossing with the first LL, there may occur other cross- 
ings with other LL. For these cases the determination of 
the extension of the non-homogeneous phase should be 
carried within a more precise method. 

The values of the transition density pl ross , and respec- 
tive proton fraction Y p and pressure P m are given for 
stellar matter under different magnetic field intensities 



in Tables II (without AMM) and III (with AMM). For 
convenience we have plotted the results without AMM in 
Fig. 4. 

For B — 0, we can see that the inclusion of the 6- 
meson reduces the transition density, and the correspond- 
ing proton fraction and pressure for NLW models, while 
for DDRH models the transition density increases and 
the corresponding proton fraction and pressure decrease. 
For DDRH models the crossing is occurring at three times 
larger proton fractions and twice the pressure. 

For a finite magnetic field it was shown in [12] for the 
models TW and TM1 (a parametrization similar to NLp) 
that the transition density increases with the magnetic 
field. We confirm this result when the S meson is included 
and summarize the conclusions: 

a) the proton fraction at the crossing increases mono- 
tonically with the magnetic field intensity, and for B* = 
5 x 10 5 it reaches the values 0.36 and 0.4 respectively for 
NLW and DDRH models. The proton fraction does not 
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FIG. 4: (Color online) Physical quantities at the inner edge 
of a compact star in the presence of strong magnetic fields: 
a) density, b) proton fraction and c) pressure. 



depend on the inclusion of the 5 meson and NLW models 
predict systematically a smaller proton fraction; 
b) the transition density is also increasing with the in- 
crease for B* < 5 x 10 4 . For smaller fields a decrease 
on the transition density may occur as it is seen with 
the DDRH models. The smaller the field, the larger the 
number of bands, defined by the filling of a new LL, will 
be crossed by the /3-equilibrium EOS. As referred before, 
this situation requires a more carefully study. For the 
smaller fields the transition densities for NLp and NL<5 
are almost coincident and only above B* = 2 x 10 5 does 
the NL<5 predict a larger transition density. Looking at 
the Fig. 1 it is seen that the symmetry energy for these 
two models start to differ at the saturation density and 
this may explain the transition density behaviour. How- 
ever the transition pressure increases faster for NL<5, and 
for B* = 5 x 10 4 NLp and NLS have the same transition 
pressure. For DDRH models, DDRLW predicts smaller 
densities for B* < 3.5 x 10 5 . The change of behaviour 
for the larger fields seems to be related with the change of 
curvature of the symmetry energy at ~ 0.12 fm~ 3 , above 
which the symmetry energy for DDRH<5 becomes much 



stiffer. It is interesting to notice that for these models 
only above B* = 5 x 10 4 do the transition pressure and 
density increase with respect to the B=0 values; 
c) even though the proton fractions and transition densi- 
ties may be similar, all models predict different pressures 
at the transition. NLW models have larger pressures for 
B* < 3 x 10 5 , NL<5 being the one with the larger values. 
For the DDRH models and below B* < 3 x 10 5 , it is the 
model including S which has the smallest pressure, about 
half the value calculated with NL5. 

For B* = 5 x 10 5 , the transition density is above the 
saturation density, p ~ 0.16 — 0.17 fm~ 3 for NLW models 
and to p ~ 0.19 - 0.23 fm~ 3 for DDRH models. For this 
high field both models with the S meson predict larger 
transition densities, pressures and proton fractions due to 
the stiffness of the symmetry energy at those densities. 

TABLE III: Same as Table II, but with AMM included. 



B* 


Model 


pr ss (fm- 3 ) 


Y P 


P m (MeVfnr 3 ) 


5 x 10 4 


NLp 


0.079 


0.079 


0.570 




NL<5 


0.079 


0.075 


0.558 




TW 


0.078 


0.097 


0.606 




DDRH<5 


0.080 


0.098 


0.498 


10 5 


NLp 


0.093 


0.152 


1.274 




NL<5 


0.094 


0.149 


1.329 




TW 


0.091 


0.178 


1.281 




DDRHW 


0.088 


0.183 


1.108 


2 x 10 5 


NLp 


0.091 


0.153 


1.206 




NL<5 


0.108 


0.260 


1.819 




TW 


0.092 


0.305 


1.355 




DDRHW 


0.088 


0.313 


1.198 


5 x 10 5 


NLp 


0.134 


0.434 


1.693 




NL<5 


0.143 


0.442 


2.245 




TW 


0.149 


0.488 


2.772 




DDRH(5 


0.170 


0.504 


4.356 



In Table III we show the same data given in Table II 
but including the AMM in the calculation. The conclu- 
sions are similar: for NLW models, the transition den- 
sity and the corresponding pressure increase when the 
(5-meson is included, whereas for DDRH models the inclu- 
sion of the isovector-scalar meson decreases the transition 
density and the corresponding pressure for B* < 5 x 10 5 . 
In summary, the S meson gives rise to a larger crust for 
NLW models and, except to very large fields, to a smaller- 
crust with DDRH models. 

The eigenvector associated with the eigenvalue of the 
free energy curvature matrix defines the direction of the 
instability and tells us how does the system separate into 
a dense liquid and a gas phase. It was shown in [13, 27] 
that in the absence of the magnetic field the direction of 
instability favors the reduction of the isospin asymmetry 
of the dense clusters of the system, and increases the 
isospin asymmetry of the gas surrounding the clusters, 
the so called distillation effect. This effect is represented 



6 



in Fig. 5 where it is seen that for the B — curve the 
fraction Sp~/5p~ is larger than p p /p n below y p = 0.5 
and the other way round above. In this figure we show, 
respectively for NLW models (left) and DDRH models 
(right), the fraction 5p~ /Sp~ as a function of y p for a 

fixed baryonic density, p = 0.06 fm -3 , chosen inside the 
instability region. 

For NLW models and for the two largest fields consid- 
ered the spinodal region contains a single Landau level 
and the curve varies smoothly starting at Sp~ /5p~ ~ 1.5 
for NLp and ~ 1.62 for NL6. We point out the very 
large value of this fraction, always above 1. The mag- 
netic field favors a strong increase of the proton fraction. 
For y p > 0.5, NLp and NL<5 behave in a similar way, 
while below this value the main difference is the larger 
8p~ /5p~ for NL<5 corresponding to a stronger distillation 
effect. 

For B* = 10 5 the spinodal has two bands, see Fig. 2, 
corresponding to the occupation of the first two Landau 
levels. The transition from one to the other is clearly 
seen with a large discontinuity of 5p~ /Sp~ at y p ~ 0.7. 
Above this y p value the curve behaves like the previous 
ones. However for y p < 0.7 the behavior is quite differ- 
ent: the curve decreases from the value at y p =0, which is 
independent of the magnitude of the magnetic field, to a 
value much smaller than the corresponding value of the 
fraction p p /p n . The same behaviour occurs for NLp and 
NL<5. The fluctuations will not drive the system out of 
the first Landau level and therefore the larger the pro- 
ton fraction, the closer the system comes to the top of the 
band and the smaller are the allowed proton fluctuations. 
For y p > 0.7 or for the larger magnetic fields the Landau 
levels are only partially filled and the fluctuations will 
never drive the system out of the corresponding Landau 
level. In summary, the effect of the 5-meson on the in- 
stability region of NLW models is to reduce the strength 
of the distillation effect of neutron-rich matter. 

For DDRH models the inclusion of the isovector scalar 
meson reduces the strength of the distillation effect, 
which is not so dramatic in these models. Including the 
AMM similar conclusions are drawn. The AMM favors 
still larger proton fluctuations because neutron polariza- 
tion stiffens the EOS. 

In conclusion, we have studied the influence of the 
isovector scalar meson on the instabilities of stellar mat- 
ter under very strong magnetic fields. The fields consid- 
ered are much stronger than the strongest field measured 
at the surface of a magnetar which is B* ~ 10 2 for SGR 
1806-20 [5]. However, the magnetic fields in the interior 
of neutrons stars could be larger and the present work 
shows how fields of the order of B* = 10 18 — 10 19 could 
affect the inner crust of a compact star. According to the 
scalar virial theorem [28] the maximum magnetic energy 
could be comparable to the gravitational energy in an 
equilibrium configuration, which would correspond to a 
value of the order of ~ 10 18 G. 

We have considered four relativistic nuclear models: 
two models with constant couplings (NLp and NL<5) and 
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FIG. 5: (Color online) 5p~ /8p~ plotted as a function of the 
proton fraction with p = 0.06fm -3 for the NLW models (NLp 
(thin lines) and NL5 (thick lines)) (left) and DDRH models 
(TW (thin lines) and DDRH5 (thick lines)) (right) and for 
several values of the magnetic fields without AMM. The frac- 
tion Pp/pn is given by the thin dotted line. 



two models with density dependent couplings (TW and 
DDRH<5). Two of the models include the <5-meson. For 
all the models, we have determined the spinodal surface 
from the curvature matrix of the free energy for different 
magnitudes of the magnetic field. It had already been 
shown [12] that the instability region could be divided in 
to several bands according to magnitude of the magnetic 
field and the number of the Landau levels occupied and 
that the presence of the magnetic field would generally 
increase the instability region. 

We have seen that the inclusion of the S meson in- 
creases (reduces) the size of the thermodynamic insta- 
bility zone for very neutron-rich matter for the NLW 
(DDRH) models when compared with the spinodal ob- 
tained without the 5-meson. These results reflect them- 
selves on the extension of the crust of a compact star 
under a strong magnetic field. 

By making a rough estimation of the transition den- 
sity at the inner crust of a compact star under a strong 
magnetic field from the crossing of the EOS with the 
thermodynamic spinodal, we have shown that the transi- 
tion density and associated pressure increases, for NLW 
models, with the inclusion of the 5 meson. On the other 
hand, for DDRH models the influence of the (5-meson is 
to decrease the transition density and associated pressure 
for fields B* < 5 x 10 5 . DDRH<5 is predicting the small- 
est transition densities, almost half the value obtained 
for NL<5. 

If we consider fields - 10 18 G (B* ~5x 10 4 ), as in- 
dicated by the scalar virial theorem, we may take the 
following conclusions: a) for conventional pulsars TW 
would predict a larger transition pressure and, there- 
fore, larger fractional moment of inertia in the crust 
than all the other models, about twice as large; b) for 
B* - 5 x 10 4 , and taking into account AMM, TW would 
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predict just a small increase of the fractional moment of 
inertia in the crust, ~ 15%, the NLW models predict an 
increase of 120% and 220%, respectively, with and with- 
out (5-meson, while DDRH<5 would predict an increase 
of 80%. If the AMM are not considered, within TW 
the fractional moment of inertia in the crust would be 
smaller than for conventional pulsars, while for DDRfM 
there would be a slight increase of 20% and for the NLW 
models an increase of ~ 100%. The different behaviours 
in magnetars and conventional pulsars predicted by the 
different models might be a possibility to impose stronger 
constraints on the EOS of nuclear matter. 

We have also investigated the direction of instability. 
If the first Landau level is only partially occupied the 
density fluctuations are such that the system evolves for 
a state with dense clusters very proton rich immersed in 
a proton poor gas. A larger proton fraction is favored 



energetically due to the degeneracy of the Landau lev- 
els. The S meson will only reduce slightly the strength 
of this distillation effect for DDRH models and increase 
it for NLW models. For particles occupying an almost 
complete Landau level, proton fluctuations arc smaller 
or forbidden and an anti-distillation effect results with a 
decrease of the proton fraction of the dense clusters. 
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